function [ out ] = fun_sigma_est(xx,theta,combos_all_unique,m_regr)

out = prod(BernsteinElement_2(combos_all_unique-1,m_regr,xx'),3)'*theta;

end